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The paper is devoted to the problem of resonances in one-dimensional disordered systems. Some 
of the previous results are reviewed and a number of new ones is presented. These results pertain to 
different models (continuous as well as lattice) and various regimes of disorder and coupling strength. 
In particular, a close connection between resonances and the Wigner delay time is pointed out and 
used to obtain information on the resonance statistics. 
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I. INTRODUCTION 

The problem of resonances, also referred to as metastable or quasi-stationary states goes back to the early days 
of quantum mechanics 0. A simple example of resonances [§| is provided by a potential depicted in Fig. [T] There 
is a wall, V — oo, for x > L and a potential V (x) — uS (x). For u — > oo, a particle of mass m has bound states at 

energies ^ {^jr) 2 l n ~ 1,2,...]. For any finite u the spectrum becomes continuous. However, the strictly stationary 
states which existed at u = oo do leave a trace in the continuum and turn into resonances. They correspond to poles 
of the scattering matrix S (E) on the unphysical sheet of the complex energy plane An alternative, more direct 

approach to the problem of resonances amounts to solving the stationary Schrodinger equation with the boundary 
condition of an outgoing wave only [H,0. Thus for the potential in Fig. [T]one has to solve the equation 

cPip , . , ~ 2 . ( 2mu\ 
—z- + ao \x\ w = k w \ a = —z— 



■^+a5(x)ip = k 2 ijj \a=—j (1) 

with the boundary condition ip (x = V) = and the outgoing wave condition tp (x) — e~' lkx for x < 0. The latter 
condition makes the problem non-hermitian: the eigenvalues for k, and for the corresponding "energies" E — h 2 k 2 /2m, 
will be generally complex. 
The solutions of Eq. ([T]) is 

,. , ( Asmk(x-L), 0<x<L 

^ J W = 1 -ikx ^ n 2 

I e lKX , x < 

Matching the function and its derivative at x = results in 

^1 — i"/kl?j tan kL = — ^ykL, (3) 

where 7 = (aL) 1 « 1. For 7 = one recovers the bound states, k n L = irn. For small 7 the solutions of (|3|) are 
obtained by iteration: 

k n L = irn (l - 7 + 7 2 ) - i {imf ~/ 2 + O (j 3 ) . (4) 

One can immediately write down the "eigenenergies ", E n = h 2 kl/2m = E„ - |r n . The real part, E n , gives the 
position of the resonance on the energy axis, whereas T n determines the resonance width. For n not too large, 
namely, n <C 7 _1 , the resonances are sharp, i.e., their width is much smaller than their spacing on the energy axis. 
This simple example demonstrates how true bound states in a closed system (u — > 00) turn into resonances, when the 
system is opened to the outside world (finite u. i.e. non-zero coupling constant 7). 

Open quantum systems can be described in terms of an effective, non-Hermitian Hamiltonian whose complex 
eigenvalues give the position of the resonances in the complex energy plane (in addition, there might be real eigenvalues 
which correspond to the bound states) . Such non-Hermitian Hamiltonians have been used for a long time in scattering 
theory, including scattering in disordered and chaotic systems pal . ,M , M l , B ■ There is a considerable amount of work 
on resonances m disordered potentials @ , [13 , , [ll , El , H , [IlTlJ El • An example of one-dimensional random 
potential is depicted in Fig. [5J The potential V (x) is zero for x < and it is infinite for x > L. In the interval 
< x < L, V (x) is a random function of x, with zero mean and some well defined statistical properties. There is 
also a barrier uS(x) at x — which allows to tune the coupling strength to the external world. For u — > 00 (closed 
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Figure 1: An illustration to the problem of resonances in a potential comprised of a 8- function barrier and a hard wall. 



system) all states are localized within the system. Two such localized wave function are schematically shown in the 
figure: tpE (x) is a state of positive energy, localized far away from the boundary x = 0, i.e. its localization center Xq 
is much larger than the localization length f. The function tpE< (x) corresponds to a negative energy state, which is 
localized essentially in a single deep potential well. When u is made finite the localized state ipE (x) will turn into 
a narrow resonance, with a width T proportional to exp (— 2xo/£), while the state ipE' (x) will remain a true bound 
state. A theory of resonances in disordered chains should consider the statistical ensemble of all possible realizations 
of V (x) and produce the probability distribution P (V). 

A quantity closely related to the resonance width is the Wigner delay time [l8j],[l9[ which is a measure of the time 
spent by the particle in the scattering region and is defined as the energy derivative of the scattering phase shift. For 
the single-channel scattering, as presented in the setup (b) in Fig. [2J the solution of the scattering problem amounts 
to finding the phase 9 (E) of the reflected wave, e - lkx + iS ( E ) ^ due to the incident wave e lkx . The corresponding Wigner 
delay time is defined as 

For a disordered system, 9 (E) and r (E) are random quantities, characterized by the joint distribution Pe,l {9, t) 
over the ensemble of realizations. There exists a lar ge b ody of work on the statistics of delay times for the scattering 
on disordered and chaotic systems [20], [21], [22], [23], [24], [25], [26]. In the presence of a sharp, well isolated resonance 
E n = E n — |r n , delay time at the energy E close to E n is given approximately by pl| 

T(E)^h ^ , (6) 

{E n -Ef + Tl/^ 

which demonstrates the intimate relation between the resonance width and the delay time. Below (section IV) we 
obtain a relation between the delay time and the resonance width distributions which is exact in the limit of weak 
coupling to the lead, and which enables us to obtain information about resonances based on the existing knowledge 
of the delay time statistics. 

Statistics of resonances and of delay times (or the closely related "dwell times") are of great interest in the physics 
of disordered media. For instance, in a disordered conductor the current carriers can be trapped for a long time, which 
lead to long tails in the decay of an electric current [27[ ■ Although our discussion is limited to "matter waves", obeying 
the Schrodinger equation, similar phenomena occur for electromagnetic waves as well. When a wave is injected into a 
random dielectric medium, it can spend there a very long time, before escaping from the sample. This phenomenon 
of long delay times has been extensively studied in experiments [28]. Resonances and long escape times might be also 
relevant to the phenomenon of "random lasing", when an active random dielectric medium without any prefabricated 
cavities, exhibits lasing above some excitation threshold (29| . 

The organization of the paper is as follows. In section [II] we introduce a tight binding model and, following [l4j).fl5| 
derive the effective non-Hermitian Hamiltonian for the resonance problem. Section IIIII is devoted to the case when 
coupling between the disordered chain and the external lead is weak. Both the tight binding model and a random 
continuous potential are treated. In section IIVI a relation between the distributions of resonances and delay times 
is obtained, and used for studying properties of the resonances under various conditions (weak and strong coupling, 
finite and infinite chain). Section fVl specializes to the case of strong disorder, for the tight binding model, using the 
locator expansion technique. 
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Figure 2: Schematic illustration of a ID disordered system with one end (x = L) closed, and the other one (x = 0) coupled to 
the lead through the (^-function barrier. The resonance problem corresponds to (a): the outgoing wave condition is imposed 
and complex values of k are found. The standard scattering problem is described in the set-up (b), where a particle with energy 
E = h 2 k 2 /2m is impinging on the system. 



II. A TIGHT BINDING MODEL FOR RESONANCES AND ITS EFFECTIVE HAMILTONIAN 

Along with the model of a continuous random potential, described in the Introduction, we also consider the tight 
binding (Anderson) model (TBM) depicted in Fig. [3] Black dotes, labeled by n — 1,2,. ..TV, designate sites in the 
chain. Each site is assigned a site energy e n chosen from some distribution q(e). The energies on different sites are 
independent of each other. Open circles, labeled by n = 0, — 1, —2, . . ., represent the perfect semi-infinite lead to which 
the chain is coupled. The lead simulates the free space outside the chain. All nearest neighbor sites of the chain are 
coupled to each other by a hopping amplitude t, and the same is true for all nearest neighbor sites of the lead. The 
only exception to this rule is the pair n = (0, 1) which provides coupling between the chain and the lead. The hopping 
amplitude for this pair is taken to be equal t' . This allows us to tune the coupling from t' = (closed chain) to t' = t 
(perfect coupling). The Schrodinger equation for the entire system (chain + lead) is a set of coupled equations: 

-ttpn+l ~ tlpn-l = Etpn {n < 0) (7) 

-tip-! - t'tpi = Expo (n = 0) (8) 

-tfh - t'lPo + ei^i = Etpx (n = 1) (9) 

-tip n+1 - ti\> n - X + e„^„ = Ei/) n (1 < n < N). (10) 

with the Dirichlet boundary condition ^jv+i = 0. Eqs. (fT|)- (fTt))) are to be solved subjected to the boundary condition 
of an outgoing wave in the lead, i.e. iji n oc exp(— ikn), for n < 0, with Refc > (the wave propagates from right to 
left). The complex wave vector k is related to E by E = —2tcosk. The complex solutions E a = E a — |r a of Eqs. 
(T7|- (fTU| yield the width of the resonances, as well as their position along the energy axes E. 

As has been explained in the Introduction, the condition of an outgoing wave makes the problem a non-Hermitian 
one. In particular, for the tight binding model (Fig. [3]) one can derive an explicit expression for an effective non- 
Hcrmitian Hamiltonian whose eigenvalues correspond to the resonances, in addition to the possible bound states. 
Using the plane wave solution ip n oc exp(— ikn) in the lead (n < 1), it is straightforward to eliminate from Eqs. 
([7|)- (fT0f all ijj n 's with n < 1 (for details see |6j), thus reducing the problem to a system of equations for the amplitudes 
ip n on the sites of the disordered chain alone (n = 1,2, ...N): 

- tip n+ i - tip n _i + e n ip n = Eip n (n = 1,2,..., N) (11) 

with the boundary conditions tJjq — ipN+i — 0. Here e n = e n for n = 2,3,...., but not for n = 1. This end site is 
assigned a complex energy 

e 1 = e 1 -tne ii , (12) 

where the parameter n = (t' /t) 2 describes the coupling strength to the outside world. Thus, the effective non- 
Hermitian Hamiltonian H . defined in (TTT|) . differs from the Hermitian Hamiltonian, H, of the corresponding closed 
system (i.e., with 77 = 0) only by the complex correction to the energy of the first site (the only site coupled directly 
to the lead), i.e., 

H = H — tr ie tk P, (13) 
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Figure 3: Resonance problem for one dimensional TBM. 



where P is the projection on site n = 1. Note that the effective Hamiltonian H depends, via k, on E. Therefore 
Eq. (fTTj) does not constitute a standard eigenvalue problem and the eigenvalues of H have to be determined self- 



consistently. We denote the complex variable E 



E 



\T. It is shown in [If 



that the resonances z Q , in the 



complex z-plane, correspond to the roots (with Imz Q < 0) of the equation 



z-a - S^z) + V te fHz) = 0, 



(14) 



where S\(z) is the self-energy for site n = 1 and k(z) is related to z by z — — 2icosfc (in [15] the variable z was 
measured in units of t). 



III. TREATING THE COUPLING TERM IN H AS PERTURBATION 

When the coupling to the lead is weak (7/ <C 1), the resonances can be obtained as small corrections to the eigenvalues 
of the closed system. For the tight binding effective Hamiltonian H, Eq. (fT5|) . first order perturbation theory with 
respect to the coupling term ~t-qe %k P gives 

E a = e a - rfti>l (1) e ik « = E a - % -T a , (15) 

where e a is the energy of the unperturbed eigenstate ip a [the former notation tp n has been changed into ip a (n) where 
subscript a labels the eigenstates] , related to k a by 



-2tcosk a ^e ik " =-g+*yi-ik. (16) 



Thus, the resonance width is 



T a = 2rpfc (l)sinfc Q = 2r^ (1) ^t 2 - \e%. (17) 

In addition to the imaginary correction, — ?X/2, there is also a real-valued correction, i.e. E a — e a + rjip 2 , (1) e a /2. 
This small energy shift on the real axis is of no interest. Note that the resonances exist only for \e a \ < 2t, i.e. within 
the band of the lead. For energies outside the band only bound states exist (real E a ). 

An expression analogous to (|17[) is obtained also for the continuous case depicted in Fig. either as a continuum 
limit of (|17p or by direct application of the perturbation theory. In the latter approach, matching the internal solution 
to the outgoing wave in the lead, one obtains 



0+,/c 



tp 1 (a 



= 0+,fc) 



' (18) 



2mu 



Here ip (x,kj is the solution in the interval < x < L for the energy E — A^-, which satisfies the closed-end 
boundary condition tp — L,k\ =0 plus the condition of the outgoing wave for x < 0. For u — > 00, Eq. p^|) gives 
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the spectrum and the eigenstates of the closed system, satisfying zero boundary condition i> a (0,k a ) — 0. For weak 
coupling to the lead, g = jtt- ^> f , perturbative expansion of the above secular equation in powers of g~ l , yields 



h 2 k 2 



2m 



2mu 



4>k (0,k a )' 



where ip^ (x, k) 



Ok 



ip(x,k). Then, employing the identity (see, e.g., Refs. [25], [19]) 



1™ J_ I J~ JP V dxdE 



2m dx \ dx dE 



(19) 



(20) 



the resonance width is expressed as 



r„ = 



2m 



ft hot 



mu 



' iC (O.fcq)l 2 

2kl ' 



where ip a is the normalized to unity eigenfunction of the closed system with the eigenenergy E a — 



2m r 



(21) 
This 

expression is consistent with the exact effective Hamiltonian for the continuous open systems derived in Refs. [30l|.[3ll|. 

Thus, both in the continuum and in TBM, for weak coupling to the lead there is one-to-one correspondence 
between the resonances and the eigenstates of the closed system, and the resonance width is related to the tail of the 
corresponding eigenstate at the boundary. 

Certain simplifications occur in the limit of a semi-infinite chain. The TBM in the TV — > oo limit has been studied 
in (TBI, where the small-r asymptotics for the density of resonances (DOR) has been rigorously derived in the weak 
coupling limit (rj <C 1). DOR in the E, T-plane, for a given realization of the disorder, is given by 



p(E, r)=J2 S(E - E a )6(T - r a ), 



(22) 



where z a — E a — |T a are solutions of (fTT| . For any T ^ 0, this expression for DOR has a well defined N —> oo limit 
and no division of the sum by N is necessary, - in contrast to the usual case of the density of states (on the real axis) 
for a Hermitian problem. Note that the probability distribution of resonance width P(T) (for some fixed E) does 
not have a well defined N — > oo limit and it approaches S(T). (Indeed, for a semi- infinite chain an eigenstate will be 
localized, with probability 1, at an infinite distance from the open end and, thus, will be ignorant about the coupling 
to the external world.) Thus, the appropriate quantity to look at for a semi-infinite chain is the DOR, rather than 
the probability distribution of resonance width. This subtle point is discussed in some detail in (T^ . 

Although the general considerations in [l5j pertain to any coupling strength 77, specific results for the average DOR 
where obtained only in the weak coupling limit, where the width of all resonances becomes proportional to rj. The 
small-r asymptotics for the average DOR (p(E, T)) is [lj| (wherein the result is written in terms of some rescaled 
variables): 



(P(E,T)) 



v{E)£,{E) 
2T 



(23) 



where v(E) and £(-E) are, respectively, the usual density of states (on the real energy axis) and the localization length 
for an infinite disordered chain. Angular brackets denote averaging over the ensemble of all random realizations. This 
asymptotic (1/r)- behavior is universal, in the sense that it holds for any degree of disorder and for any — 2t < E < 2t. 

The 1 /r-asymptotics can be understood with the help of a simple intuitive argument which, in somewhat different 
versions, has appeared in [3,[l(|j[Ili],[12],[25]. The essence of the argument is that narrow resonances stem from states 
localized far away from the open boundary, say, at distance x. Such states will have an exponentially small tail at the 
boundary, proportional to e~ x /^, and the corresponding resonances will be exponentially narrow, T - e - 2x ^. The 
1/r- behavior then immediately follows from the assumption that the localization centers, x, are uniformly distributed 
in space. 

One should keep in mind that, for a long but finite chain of N sites, the (l/r)-tail will be cut off at very small T 
of the order of exp ( — 2Ay£). The extremely narrow resonances with T exp ( — 2A^/£) originate from states localized 
in the vicinity of the closed-end site n — N and they should be treated separately (see below). 



IV. RELATION BETWEEN DISTRIBUTIONS OF RESONANCES AND DELAY TIMES 



The rigorous asymptotic result of the previous section, Eq. (|23p . was obtained for a semi- infinite chain weakly 
coupled to an external lead. Things get more complicated if these restrictions are relaxed. In particular, the simple 
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relation between the resonance width and the behavior of the corresponding eigenstate of the closed system [Eqs. (fTT)) 
and (|2"Tj) ] breaks down when the coupling between the system and the lead becomes strong. In this section we discuss 
systems of finite size L and beyond weak coupling limit. 

In the Introduction we have mentioned the problem of the delay time r (E, L) and the corresponding phase shift 
6(E,L), for a particle of energy E impinging on a random chain of length L. We designate by Pe,l (#, t ) the joint 
probability distribution of 9 and r for perfect coupling to the lead (r) = 1, or u = 0) and relate this distribution to 
the average DOR (p (E, V)). Such relation is useful because it enables us to "transfer" the existing knowledge of the 
time delay in disordered chains |2(J,|2j|,[2j^,[23],[24],[25j into the field of resonances. To this end we introduce the 
quantity 



J. , ■■. 



< (0) 



(24) 



where ip a (x) is a normalized eigenfunction of the closed system satisfying the boundary conditions if) (0) = if) (L) 
The average density of points {E a , $ Q } in the (E, $)-plane is 



0. 



(25) 



Although (p (E, $)) is defined in terms of eigenvalues and eigenfunctions of the closed system, it can be related to the 
distribution Pe,l (9, t) which describes scattering properties of the corresponding open system. The relation stems 
from the fact that for 8 (E) = it the scattering wave function vanishes at x — 0, so that the eigenvalues E a are given 
by zeros of the function 9 (E) — tt. This observation results in the identity 



d6{E) 



dE 



S(E — E a ) 



r(E) 



(26) 



A generalization of this identity involves, in addition to the eigenvalues E ai also the eigenfunction-related quantity 
and it reads [32] 



J2HE-E a )5(®-$ C 



(27) 



where v g = h is the group velocity in the lead. This identity, upon averaging over the distribution Pe,l {@i t ) 
and using (|25|) . yields the required relation between the quantities characterizing the open and the closed system: 



(p (£,*)) 



P E ,L [e = 7T,T 



(28) 



This expression holds for arbitrary L and has a well defined L — > oo limit [cf. the discussion after Eq. (|22l) ]. Let us 
note that Eq. (|2"8"|) constitutes the strictly one-dimensional counterpart of the similar relations derived in Ref. (33| 
for the one-channel scattering from a higher-dimensional system. The results in Ref. [33j were obtained within the 
nonlinear sigma model and, thus, do not include the strictly ID case discussed here. 

Equations (|24| . ([28|) correspond to the continuous model. A completely similar treatment for the TBM yields 
precisely the same relation (|28|) [with h = 1], but with $ Q redefined as 



2 sin 2 k ' 



(29) 



and the group velocity in the lead given by v g = 4j§ = 2isinfc. 

The relation (f2"8"| is rather general. It holds for an arbitrary L, for any degree of disorder, and it is applicable to 
lattice models as well as to continuous ones. However, to employ this relation for the resonance statistics problem 
one more step is needed, namely, a relation between $ a and r a . For the weak coupling case such a relation has been 
derived in the previous section for both the TBM [Eq. ([TT)) ] and the continuous potential ([2"Tjh The two expressions 
can be unified into a single formula 



T 



hv a 



(T«l) 



(30) 
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where T is the transmission coefficient through the potential barrier separating the lead from the chain. The latter is 
realized by a ^-function potential in the continuum or by the weak hopping link t' in the TBM, as described previously, 
so that 



T 



continuous model 



l + (mu/h 2 kf ' 

n \Y!f k 2lr , TBM 

(1— i)) sin 2 k ' 



Note that the linear relation (|50|) between T a and $ Q is valid only if T is small (weak coupling). With the help of 
(|30p one can map the density (p (E, $)) in the (E, $)-plane, Eq. (|28p. onto the average DOR in the (E, T)-plane: 

HT 2 ( h\ 
( p (E 7 r)) = —P E , L (6 = n,T = T-). (32) 



p3 ^'^ \ ' p 

This formula relates the average DOR to the delay time statistics. For a weak Gaussian white noise disorder and 
L £ the distribution Pe,l (#> t) does not depend on and has the following form [24J: 

PeMS , t) . _^ e -W, + _L e -W 2 ^^ e -,( 1+ .-)«_^_ stoh « H , w , /2 (2) , (33) 

where To = £/u fl an d ^i,is/2 is the Whittaker function (the same result is obtained for the weak correlated disorder 
[12] )■ Expression (|3"3"|) , via (132p . immediately yields the corresponding DOR. 
In the limit L/£ — > oo, r fixed, Eq. (f3"3")) reduces to 



Pe,»(^)4V T . ( 34 ) 

Z7T T 

so that 



(p(£;,r))^^e- r / r o, (35) 

where Vq = (irhvg) 1 is the density of states in the lead per unit length and 

To = T^. (36) 

Eq. (|35p coincides with the former result (|23[) for T/ro <C 1 and, in addition, gives an exponential suppression of the 
resonance density for T/Tq > 1 (the exact density of states v in (|23|) reduces to v$ in the weak disorder limit). 

For finite size chain (L/£ « 5) and weak coupling (T .0004) the distribution P(logT) = (T/vL) (p (E, T)), 
calculated from (|32 |) . ((33|) is presented in Fig. IIVI (solid line). For comparison, a numerical Monte-Carlo simulation 
was performed for the TBM at the energy close to the unperturbed band edge (dots in Fig. HV[) . The agreement is quite 
good. Let us discuss the example in Fig. IIVI in more detail. First, the exponential factor exp [— r/To], which suppresses 
the large-r probability, is present in both the semi-infinite, Eq. (|35|) . and finite-i, Eq. (|33|) . case. For T's smaller 
than the characteristic value To, one can distinguish two regimes. In the intermediate regime, e~ 2L ^ <C T/Tq < 1, 
the behavior P(logT) m const [i.e. (p(E,T)) ~ 1/r, Eq. ([33)1 ] is valid, since the opposite closed boundary of the 
system has not yet come into play. On the contrary, the regime of very narrow resonances, T <C Toe~ 2L ^ , is strongly 
affected by the boundary x — L. These resonances are associated with the eigenstates localized close to this boundary 
and are described by the nearly log-normal tail of the distribution. 

Although the regime of the narrow resonances, V <C roe _2L /^ is contained in the analytical expressions (|3"2")) . (|33[) . 
it is worthwhile to give an independent, more direct derivation. Let us recall that the localization length £ is defined 
for an infinite system and, in this limit, it is a self-averaging quantity. In a long but finite size chain (L ^> £) the 
localizaton length, or more precisely its inverse (the Lyapunov exponent A) is a fluctuating quantity with nearly a 
Gaussian distribution (see e.g. (25| and references therein) 

P x{ X;L) = Me- L ^- 1 ) 2 / 2 . (37) 
V Ztt 

The tail of the extremely small T's is related to the eigenstates localized near the closed end of the system, x = L, for 
which $ Q fa $ e _2iA , where the pre-factor <&o ~ is of minor importance. Then, neglecting the pre-exponential 
factor, the probability for $ Q -C <i>oe~ 2L /£ decays like 



P (<&) oc exp 



(ln($/$ ) + 2L/Q 2 



(38) 




Figure 4: Weak coupling to the lead. The numerical simulation (dots) was done for the TBM (f7))- (|10[) with i = 1, E = —1.9, 
r\ — 0.001 and system length N = 401 sites. The uncorrelated disorder with a box distribution was implemented, which 
produced localization length £ = 77.6 sites. The analytical curve was calculated according to Eqs. (|32[) and (|33l) . 



Using (|30| . one obtains the log- normal cutoff of the DOR 



(p(E,T)) cxexp 



(ln(r/r ) + 2£/Q s 



« e- 2i /« 



(39) 



Similar cutoffs for the delay time and the average DOR have been derived in Refs. [25| and [ll| respectively. 

So far the discussion was limited to the weak coupling case, when a simple relation between $ Q and r a [Eq. (|30p ] 
could be rigorously derived. When the coupling parameter T increases and approaches unity, the relation (|30[) ceases 
to be quantitatively accurate and turns into an order of magnitude estimate T a ~ hv g <fr a . This relation is physically 
reasonable for narrow, isolated resonances. Such resonances stem from the eigenstates (of the closed system), which 
are localized far away from the open boundary x — 0, and their width is much smaller than the mean level spacing. 
One can then trace a particular resonance, i.e. its width T as a function of the increasing coupling strength T, without 
worrying about other resonances. It is therefore intuitively clear that the small- T result, Eq. (|30p , can be qualitatively 
extrapolated up to the perfect coupling limit T = 1. 

One can support the above argument by a more elaborated analysis. Consider the formal solution of the one-channel 
scattering problem at energy E close to a narrow isolated resonance z r — E r — ^T. For small E — E r , using general 
analytical propertiesjaf the scattering amplitude in the complex energy plane, the solution in the lead can be expanded 
as (see, e.g., Ref. 



i/f(x,E) = a[E-E r + ^r\e ik * + a* (e — E r — ~T ) e~ ikx , x < 0, 



(40) 



where, by identity (|20p , the complex constant a satisfies (up to small corrections) 



1 



\tp (x)\ dx. 



Energy E a at which ip (x — 0, E a ) = is the eigenenergy of the closed system, and by (fH))) 

E^ — E r — r. 

2 Re a 

For the corresponding eigenfunction of the closed system, using P0|) , (I41|) in the definition 

2 



one obtains 



2hv 



Ima 
Re a 



1 



> 



r 



2hv„ 



(41) 



(42) 



(43) 



Both and (|43p are meaningful as long as [E a — E r ) < T, i.e., (jrjp ) 2 < 1, since otherwise the linear expansion P0|) 
is not valid and higher orders should be included. With this reservation, Eq. (|43|) relates narrow isolated resonances to 
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Figure 5: Resonance width distribution for the perfect coupling T = 1. The analytical curve is calculcated according to Eqs. 
(I32ll , (|33|) with To given by (145|> . The numerical simulation was done as explained in Fig. IIVI with the TBM parameters i = 1, 
E — —1.9, rj = 1 and the uncorrelated disorder resulting in the localization length £ = 77.6 sites. The system length is (a) 
N = 401 sites and (b) N = 601 sites. The scatter of the numerical data for the extremely small values of F/Fo occurs because 
of the insufficient numerical statistics. 



the well localized eigenstates of the closed system. However, contrary to (|50"]) , relation (|4U|) is not deterministic, since 



it depends on the phase of a (which is random for weak disorder). Replacing the unknown coefficient 
by a phenomenological constant leads to 

T a = 2phv g $ a . 

With the relation (|44|) at hand, all the steps done for the weak coupling can be repeated, and Eqs. 
(155)1 apply with the transmission coefficient T replaced by 4/3 and the characteristic value T [Eq. (]3"6"]) ] 



Ima\ 2 
a / 



r„ 



VgH 



4/3^ 



>Rc 



(44) 

EH),® and 
redefined as 

(45) 



In the present case, however, the DOR obtained from Eqs. (|32|) . (|33|) is valid only for T <C To, since otherwise the 
isolated resonance approximation implied in the above argument is not applicable. 

The above approximation was compared to the numerical simulation for the perfect coupling to the lead, Fig. [5] 
In both cases shown in Fig. L/£ « 5 and L/£ w 8, the same fitting value f3 « 0.68 was used. As expected, a good 
agreement between the numerical simulation (dots) and the analytical result (solid line) is obtained only for T < To 
(the deviation for the extremely small T's is due to the numerical under-sampling). 



V. STRONG DISORDER 



In this section we consider the case of strong disorder, when the hopping amplitude t is much smaller than the 
characteristic width W of the site energy distribution q(e). For a semi-infinite chain the problem was considered in 
[l4| . making use of a recursion relation for the self-energy. Here we employ the locator expansion, i.e. perturbation 
theory in i, which is the appropriate tool for strong disorder (13, [36]. Our treatment is not restricted to a semi- 
infinite chain and, in particular, we address the question of the cutoff of the (l/r)-tail in a chain of large but finite 
N. Furthermore, no restriction on the coupling strength rj is imposed in our treatment. 

For t = the Hamiltonian ([IT?]) corresponds to uncoupled sites and its eigenvalues coincide with the site energies €j 
(j = 1, 2, ....N). When t is switched on, some of these "unperturbed" eigenvalues acquire a complex correction, due 
to the last term in (|13p . and thus describe resonances. Our purpose is to find the imaginary part of this correction, 
in the leading order in t. (The small correction to the real part, e^, introduces an unessential shift on the real axis of 
the complex energy plane and will be ignored). We designate the complex energy E by z = E — and look for the 
solutions, Zj, of Eq. (|14p . which we rewrite as 



z — l\ — S\{z) = 0, 



(46) 
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with 

h = ei -rjte 1 ^. (47) 

In order to see the mechanism by which the unperturbed solutions, Zj = tj, acquire an imaginary correction, we 
employ the locator expansion for the self-energy Si (z). It can be represented diagrammatically as a sum over all 
paths which start at site 1 and return to this site only once [35|.[36|. An example of such a path is drawn in Fig. [5] 
This path goes from site 1 to 2, proceeds from 2 to 3 and returns back to 1 . This path contributes to Si a term 
tgitg%tg%ty where g n = (z — en) -1 is the Green's function (the locator) for an isolated site n. Thus, the general rule 
is that to a line connecting a pair of sites one assigns the number t, while to a site n the corresponding locator is 
assigned. By inspecting Eq. (|46[) it becomes clear that an imaginary correction to the unperturbed solution z^ is 
produced by paths, in the Si - expansion, which connect site 1 to site j. Indeed, site j has no direct knowledge about 
the connection to the outside world: this information must be transmitted to it from site 1, via all intermediate sites. 
To leading order, it suffices to keep the shortest path. For site 3 this is the path in Fig [51 Generalization to an 
arbitrary site j is obvious and results in a path of (j — 1) loops which brings in a factor i 2 ^ -1 ). This path produces 
the imaginary part of Zj , which is calculated from Eq. (1461) : 

1 1 

Im Zj = t 2 ^- 1 ) TT - Im -. (48) 

fc=2 ^ _ Cfc ) ^ ~ 61 

Since only the leading term (in powers of t) is kept, we have replaced in all the locators z by z^ = ej. For the same 
reason, k{z) in the expression (|4"7T) can be replaced by k(E). From the relation E = — 2tcosk it follows that 



gi=ei+ ^(l- i v i -5j- (49) 

Note that the imaginary part in (|4"5)) exists only for \E\ < 2t, i.e. only bound states in this energy interval (in a closed 
chain) turn into resonances upon coupling the chain to the lead (the same energy interval has already been identified 
in Sec. IHIj) . Eigenstates beyond this energy interval remain strictly bound states. Substituting (|49|) into (|48| and, 
again, keeping only leading terms in t, one finally obtains: 



The DOR in the (E, r)-plane is given by 



p(E,T) = J2HE-e j )S(T-r j ). (51) 

i=2 

Since the small shift of the eigenvalues along the real axis is of no interest, we have set Ej = ej in Eq. (jSTj) . For a 
fixed j, resonance width Tj depends on the energies of all previous sites, k — 1, 2, ... ,j — 1, but not on ej. Therefore, 
the two (5-functions in (|51l) are statistically independent, so that upon averaging 

(p(E,T)}= q (E)^f28(T-T j )^, (52) 

where, in the strong disorder limit, the site energy distribution function q (E) coincides with the density of states per 
site in the closed system. To avoid cluttering the notation we set E = (middle of the band) and 77 = 1 (perfect 
coupling). (Extension to arbitrary E and rj requires some obvious minor modifications.) 
For this case 

( P (E = 0, r)> = q (o) J2 (r at J] I) J ■ (53) 

It is convenient to define a random variable 
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Figure 6: Calculation of the self-energy in the locator expansion. 



This is a sum of independent random variables with the average value ( In Q \ = J deq (e) In ^ = a and variance 



In 



2 £ 



tW 



< 



of — j3. For instance, for the Anderson model, when e is uniformly distributed within a window 

e < ^Ip, one has a = (in ^ — l) and j3 = 1. Note that, in the strong disorder limit, a coincides with the inverse 
localization length (the Lyapunov exponent) [37] . 

Since N is a very large number, most of terms in (|53[) correspond to large j, so that Aj has a Gaussian distribution, 
Pj (A), with the average value (Aj) = —2 (j — 1) a and variance (AA|) = 4/Jj, i.e. 



1 



exp 



(A + 2jar 
8J/3 



(55) 



where (j — 1) was replaced by j. Eq. (fSU)) then yields 



= 



,r))=g(0)X; / 8(T-2te A )P 3 (A)dA 



9(0) 



J=2 

E 

J=2 



1 



exp 



(ln(r/2t) + 2j/Q 2 
Si/3 



(56) 



where £ = l/o is the localization length in the middle of the band (E = 0). The lower limit of summation, j = 2, 
should not be taken literally and it is of no importance, since for small resonance width T the sum is dominated by 
large- j terms. 

For narrow (but not too narrow) resonances, when 1 <C — In (T/2t) <C N/£, the sum is dominated by terms with j 
near jo ~ — (£/2) In (T/2t) 3> 1. Then, the sum in (|S"6")) can be approximated by an integral and 



(p(E = 0,T)) 



g(0)g(0) 

2r 



(57) 



in agreement with the universal result in Eq. (|23p . This 1/r behavior is cut off sharply for very narrow resonances, 
such that — ln(r/2t) ^> N/£. These resonances stem from states which are localized in the vicinity of the sample 
boundary at j — N. The sum (|5"6")) is then dominated by the last term, i.e. 



(p(E = 0,T)) 



9(0) 



r 



exp 



(In (T/2f) + 2N/CY 
87V/3 



(58) 



i.e. for — ln(r/2f) > 2N/£ > the DOR rapidly (faster than any power of T) approaches zero with decreasing T. This 
kind of log-normal tails are well known in the theory of disordered electronic systems [27]. 

It is instructive to compare the strong disorder result, Eq. (|58p. with the expression (|31?)) which was derived in 
the opposite case of weak disorder. The main difference between the two expressions, besides the fact that in (|31?)) 
the pre-exponential factor has not been written down, is that the exponent in (|39p contains the single parameter 
L/£, whereas (|58|) depends in addition on the parameter (3/a [indeed, 8Nf3 can be written as 8(N/£) (/3/a)]. The 
parameter (3/a is a non- universal number which depends, for instance, on the chosen distribution for the site energies, 
q (e). The same situation is well known to occur in the study of the transmission coefficient T through a disordered 
chain of length L. The distribution of InT is Gaussian. If the disorder is weak, then there is a universal relation 
between the mean and the variance of InT (single parameter scaling). On the other hand, for strong disorder the two 
become independent of one another (two parameter scaling) (38[. 



VI. CONCLUSION 



Statistics of resonances in disordered one-dimensional chains is a formidable problem which does not easily lend 
itself to a rigorous analysis. In this paper we have reviewed some of the existing results and have extended them 
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in various directions. We consider both a continuous random potential and the tight binding lattice model, and we 
tackle a variety of different cases, differing by size L of the chain, by strength of the disorder or by coupling strength 
between the system and the external world. There is no efficient universal method for treating the problem in its 
full generality. Different techniques turn out to be appropriate in different regimes. In particular, we presented in 
some detail the method of locator expansion, most suitable for the strongly disordered lattice model. On the other 
hand, for weak disorder we were able to use some known rigorous results for the Wigner delay time problem to obtain 
information on resonance statistics. 

BS is indebted to H. Kunz for previous collaboration on the subject. We acknowledge useful discussions with A. 
Comtet, J. Feinberg and C. Texier. 
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